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Abstract 

We obtain an exact simple solution of the Einstein equation de- 
scribing a spherically symmetric cosmological model without the big- 
bang or any other kind of singularity. The matter content of the 
model is shear free isotropic fluid with radial heat flux and it satisfies 
the weak and strong energy conditions. It is pressure gradient com- 
bined with heat flux that prevents occurrence of singularity. So far 
all known non-singular models have non-zero shear. This is the first 
shear free non-singular model, which is also spherically symmetric. 
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Since Senovilla's discovery [1] of an exact singularity free cosmological so- 
lution of the Einstein equation representing a perfect fluid with the equation 
of state p = 3p (and subsquently in the same framework the one [2] with 
p — p), it is now being recognised that the singularity theorems [3] can not, 
as generally believed earlier, prevent occurrence of non-singular cosmological 
solutions satisfying all the energy and causality conditions. And there is no 
conflict with the theorems in this. The theorems became inapplicable because 
one of the assumptions, existence of closed trapped surface, is not respected 
by these solutions and its violation does not entail any unphysical behaviour 
for the matter content. This assumption was however always a suspect [4] 
but this fact was not fully appreciated in absence of a non-singular solution. 
The Senovilla solution did this signal service of dispelling the folklore belief. 

A large family of non-singular cosmological models [5] and its generaliza- 
tion with heat flux [2] has been considered but they are all cylindrically sym- 
metric (see an excellent recent review [6]). For practical cosmology, spherical 
symmetry is however more appropriate. It is therefore pertinent to seek 
spherically symmetric non-singular models. The first model of this kind was 
obtained by one of us [7] which has imperfect fluid with heat flux (note, the 
expression for 9 should have a negative sign before it) and it satisfies all the 
energy conditions and has no singularity of any kind. It was obtained by 
letting one of Tolman's solutions [8] expand. The solution has a free time 
function which can be chosen suitably to have non-singular behaviour for 
physical and kinematical parameters and there exist multiple such choices. 
It is also possible to have a non-singular model with null radiation flux [9]. 
These models are both inhomogeneous and anisotropic and have the typical 
behaviour, beginning with low density at t — > — oo, contracting to high den- 
sity at t — and then again expanding to low density as t — > oo. Nowhere 
any physical parameter diverges. 

In the Raychaudhuri equation [10], which governs cosmological dynam- 
ics, it is acceleration (pressure gradient) and rotation (centrifugal force) that 
counteract gravitational collapse. In cosmology, there is absence of overall 
rotation and hence for checking collapse to avoid singualrity presence of ac- 
celeration becomes necessary. All the known non-singular models [1-2,5,7,9], 
are not only accelerating but also shearing. Though shear acts in favour 
of collapse in the Raychaudhuri equation but its dynamical action through 
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tidal acceleration makes collapse incoherent which acts against concentration 
of large mass in small enough region. This would ultimately work against 
formation of compact trapped surface. 

Raychaudhuri [11] in one of his recent theorems establishes that the nec- 
essary condition for non-singular cosmological model is that space average of 
physical and kinematical parameters must vanish. That means the parame- 
ters must depend upon space variables. The space gradient of expansion is in 
vorticity-free spacetime given by space divergence of shear and heat flux [12]. 
Hence for non-singulaity presence of atleast one of them is necessary. The 
Ruiz-Senovilla family [5] of non-singular cylindrical models is the example of 
presence of shear without heat flux. It can be shown that presence of shear 
is in general essential for perfect fluid G-2 symmetric non-singular models 
[13]. The spherical non-singular model [7] has both shear and heat flux. 

Then the question arises, could heat flux alone, of course combined with 
pressure gradient, avoid singularity? This is what we wish to demonstrate in 
this letter by obtaining a simple non-singular solution which describes an in- 
homogeneous shear-free spherical model filled with isotropic fluid and radial 
heat flux. The model satisfies the weak and strong energy conditions as well 
as has a physically acceptable fall off behaviour in both r and t for physical 
and kinematic parameters. Again there is a free time function which can 
be chosen suitably to give non-singular behaviour to model and there exist 
multiple such choices. 

The metric of the model is given by 



ds 2 = (r 2 + P) 2n dt 2 - (r 2 + P) 2m [dr 2 + r 2 (d9 2 + sin^dcj) 2 )} (1) 



Here P — P(t) which can be chosen freely. The Einstein field equation for 
perfect fluid with radial heat flux reads as 



where 



2n = 2m ± V8m 2 + 8m + 1 



in particular 




(2) 
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Rik - ^Rgik = -[(j> + P)uiu k - pg ik + -(ftUfc + q k Ui)\ (3) 

where we have set 8nG/c 2 = 1, = 1 = —qiq\qiU l = 0, p,p denote fluid 
density and isotropic pressure, and g« is the radial heat flux vector. 
From eqns. (1) and (3) we obtain 



3m 2 P 2 A 3P+(m + l)r 2 
n — — A-m - - 

1 (j.2 _|_ pyn+2 (^2 _|_ pym+2 ' 



V = 



. 2 + p )2ra+2 [2(r 2 + P)P+(3m-2n-2)P 2 ] 
__l__[( m + n)P + n V], 



4m(n + l)rP 

1 ~ (j.2 _|_ V J 

and the expansion and acceleration are given by 

3mP . nr 

We have freedom to choose the function P(t) which could be chosen suit- 
ably to give non-singular behaviour to the above parameters. As a matter of 
fact there exist multiple choices, for instance P(t) = a 2 + b 2 t 2 , a 2 + e~ bt , a 2 + 
b 2 coswt, a 2 > b 2 . For all these choices it is clear that all the physical and 
kinematic parameters remain regular and finite for the entire range of vari- 
ables. Note that it also admits an interesting oscillating behaviour in time in 
which the model oscillates between two finite regular states. The first case of 
oscillating non-singular model [14] was recently considered in the spherical 
family [7]. The oscillating non-singular models are quite novel and interest- 
ing of their own accord. 



In non-oscillating case, all the parameters given above tend to zero as 
r — > oo or t — > ±oo. The universe begins with low density and contracts to 
maximum density and then again expands to low density without ever be- 
coming singular. This is a typical behaviour for non-singular models [1,2,7]. 
However in the oscillating case, model oscillates in time between two regular 
finite states, and the parameters fall off to zero as r — > oo. This is how 
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oscillatory and non-oscillatory singularity free models differ from each-other 
in their global behaviour. 

It is obvious from the simple expression for the metric that spacetime is 
causally stable. For verification of the energy contions, we will have to find 
the eigenvalues of the energy momentum tensor, which are given as follows: 

[^(p-p + D), l -(p- p -D), p, p],D 2 = (p + p) 2 -Aq 2 . (6) 

Note that in all the above expressions there would in view of eqn. (2) 
be relative dominance of the term of (r 2 + p)- 2 ( m + 1 ), The weak and strong 
energy conditions would require p > 0, D > 0, p + p + D > 0, 2p + D > 0. It 
can be easily verified that these conditions would hold good for the choices 
for P(t) given above. The dominant energy condition which would require 
p > p cannot however be satisfied as it is clearly violated for large r. Thus 
the model satisfies the weak and strong but not the dominant energy condi- 
tion. 

We have thus obatined a spherical model with isotropic pressure fluid and 
radial heat flux without the big-bang or any other kind of singularity. This is 
the first shear free non-singular model. It is inhomogeneous but isotropic. It 
is heat flux that combines with pressure gradient to avoid singularity. From 
the point of view of realistic cosmology merits of the present model are its 
isotropy and spherical symmetry. 

Apart fron the first Senovilla model [1] and a large family of cylindrical 
non-singular models [5,2], there now also exists a large family of spherical 
non-singular models [7,9] and to that the present one adds a novel family 
of shear free models. Even though it does not satisfy the dominant energy 
condition, it is a very simple and interesting model. It is remarkable to note 
that for the first time cosmic singularity has been avoided in absence of shear. 
There has been interesting cases of cosmological models [6], for instance [15], 
which do not satisfy all the energy conditions, yet deserve consideration for 
their other remarkable properties. The present model is simple, sheer free 
and isotropic and hence is interesting enough. Above all it is a very simple 
spherical model and thus also points to an important fact that non-singular 
cosmological solutions are no longer isolated but could occur more generally 
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even in spherical symmetry. 
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